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INTRODUCTION 


This  Is  the  second  "final"  report  Issued  under  this  contract.  It  con¬ 
tains  reports  In  summary  form  of  Investigations  completed  and  new  Investiga¬ 
tions  Initiated  since  the  final  report  dated  March  30,  1962.  The  aim  of  the 
program  has  remained  as  the  vindertaklng  of  theoretical  research  to  determine 
suitable  electromagnetic  models  of  the  Ionosphere  (and  earth)  Which  can  be 
used  to  predict  terrestrial  propagation  observables  at  frequencies  below 
3  kc/s.  (Extremely-Low-Frequency  band) . 

Liaison  has  been  maintained  with  experimental  VLP/BLF  propagation  piro- 
grams  and  related  theoretical  research  at  the  National  Bureau  of  Standards, 
Central  Propagation  Laboratory,  Lincoln  Laboratory  of  the  Massachusetts  Insti¬ 
tute  of  Technology,  Air  Force  Cambridge  Research  Laboratories,  Radio  Corpora¬ 
tion  of  America,  the  U.S.  Navy  Electronics  Laboratory,  U.S.  Navy  Underwater 
Sound  Laboratory  and  the  Naval  Research  Laboratory,  Stanford  Research  Insti¬ 
tute  and  Stanford  University,  University  of  Rhode  Island,  Harvard  University, 
Uhlverslty  of  Alaska,  and  the  University  of  Texas. 

Commencing  In  December  1962,  a  special  study  was  undertaken  to  use  the 
propagation  models  developed  earlier  together  with  thunderstorm  activity  maps 
to  develop  an  ELF  noise  prediction  technique.  A  schedule  for  simultaneous 
noise  measurements  at  widely  separated  places  by  DECO  Electronics  was  set  up 
by  agreement  with  Sylvanla,  DECO  and  ONR,  NRL  and  USNUSL  personnel.  Some  of 
the  early  data  thus  gathered  was  used  to  determine  the  normalization  constant 
required  In  this  technique.  Unfortunately,  funds  on  the  contract  were  ex¬ 
hausted  before  this  task  could  be  completed. 

At  the  request  of  ONR  a  brief  study  was  made  of  the  perturbation  of  an 
electromagnetic  wave  propagating  over  the  sea  caused  by  a  submerged  conduct¬ 
ing  sphere  supposed  to  represent  a  submarine.  The  results  of  this  study  were 
presented  in  a  letter  to  A.  Shostak,  ONR,  Washington,  dated  25  August  1962. 

In  June  of  1962  Drs.  H.  Raeiner  and  J.  Oalejs  attended  the  NATO  sponsored 
ELF  meeting  In  Paris,  and  presented  papers  describing  their  work  on  the  ELF 
noise  spectrum  In  the  cavity  resonance  region,  and  ELF  propagation  below  an 
exponentially  tapered  Ionosphere. 

For  two  weeks  In  July  of  1963  Dr.  Qalejs  served  as  a  faculty  member  at 
the  NATO  sponsored  ELF  Summer  Institute  in  Bad  Honiburg,  Germany,  Where  he 
headed  the  panel  on  ELF  propagation. 
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Five  sununary  status  reports  covering  the  period  18  June  1962  to 
30  September  1963  have  been  Issued.  Several  verbal  presentations  to  Navy 
personnel  have  been  given  and,  since  the  March  30,  1962  final  report,  six 
technical  reports  have  been  Issued  bringing  the  total  to  19  In  addition  to 
two  Interim  letter  reports  on  the  ELF  noise  prediction  study. 

Most  of  the  research  completed  on  this  contract  has  appeared  In  the  open 
literature.  To  date  two  letters  and  eight  full  length  papers  have  been  pub¬ 
lished.  These  are  listed  at  the  end  of  this  report. 

This  report  consists  In  the  main  of  summaries  of  the  several  Individual 
researches  carried  out.  With  few  exceptions  the  full  details  have  been  re¬ 
ported  as  Indicated  above. 

EFFECT  OF  SEA  SURFACE  MOTION  ON  SUB-SURFACE  FIELDS 

This  project,  whose  motivation  was  described  In  the  final  report  of 
30  March  1963,  vras  completed  with  the  Issue  of  a  report  which  appeared  as  a 
paper. ^  In  this  study  attention  was  focused  on  the  fluctuations  In  the  sub¬ 
surface  horizontal  electric  field  (pointing  In  the  direction  of  wave  motion) 
caused  by  surface  waves.  Extensive  numerical  computations  of  the  mean  square 
deviation  of  this  field  component  were  carried  out  for  unidirectional  surface 
waves  of  random  height  having  an  assumed  Gaussian  correlation  function.  The 
direction  of  wave  travel,  RMS  wave  height  and  correlation  length,  frequency 
of  the  source  and  depth  of  the  observer  were  varied.  With  this  many  para¬ 
meters  the  results  cannot  be  presented  succlntly,  however,  a  few  general  ob¬ 
servations  will  be  noted  here.  The  mean  square  deviation  (MSD)  of  the 
fluctuations  decrease  with  increasing  observer  depth  and  frequency,  and  for 
small  and  large  correlation  distance. 

elfAlf  propagation  below  a  tapered  anisotropic  ionosphere 

In  the  final  report  dated  30  March  1962,  mention  was  made  of  some  pos¬ 
sible  effects  of  D-reglon  Ions  on  ELF  propagation.  A  great  deal  of  analyti¬ 
cal  work  and  numerical  computation  has  been  done  since  then.  Attention  has 
been  directed  to  cases  >^ere  the  earth's  magnetic  field  Is  horizontal.  Two 


1.  Winter,  D.F.,  "Low-Frequency  Radio  Propagation  Into  a  Moderately  Rough 
Sea,"  Jour.  Rea.  Nat.  Bur.  Stds.  67D,  551,  Sept. -Oct.  1963. 
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reports  have  been  Issued  and  submitted  £or  publication.  ’  Figure  1  Is  a 
plot  of  the  dominant  waveguide  mode  attenuation  constant  for  various  direc¬ 
tions  of  propagation  In  geomagnetic  equatorial  regions.  Figures  2  and  3  show 
the  normalized  dominant  mode  propagation  constants  In  the  frequency  region 
5  to  50  cps.  Sizeable  non-reclproclty  Is  evident  between  EN  and  WE  direc¬ 
tions  of  propagation.  As  seen  from  Figure  2,  very  large  differences  In  at¬ 
tenuation  rates  are  predicted  for  daytime  propagation  at  frequencies  from 
5  to  15  cps.  At  the  present  time  there  Is  no  direct  experimental  data  avail¬ 
able  to  check  these  predictions  against.  This  predicted  non-reclproclty  may 
be  the  reason  for  the  observed  diurnal  variations  In  the  frequencies  of  the 
peaks  In  the  ELF  noise  spectrum  observed  by  Balser  and  Wagner  and  others.  No 
detailed  Investigation  has  been  made  Into  this  possibility.  However,  there 
Is  some  reported  data  on  non-reclproclty  In  WE  and  EW  propagation  at  VLF. 
Figure  4  shows  a  plot  of  the  predicted  daytime  rates  from  reference  3.  The 
measured  data  are  In  general  agreement  with  the  predictions  but  exhibit  about 
a  half  db/1000  km  greater  attenuation. 

An  Investigation  Into  the  correctness  of  using  a  plane  wave  surface  Im¬ 
pedance  In  the  determination  of  the  mode  constants  at  VLF  was  Initiated 
towards  the  end  of  this  reporting  period.  It  was  hoped  that  this  would  shed 
some  light  on  the  problem  of  choosing  a  height  for  the  lower  edge  of  the 
Ionosphere . 

PREDICTION  OF  ELF  NOISE  LEVEL 

The  object  of  this  study  Is  to  predict  the  diurnal  and  seasonal  RMS 
noise  level  for  any  point  on  the  earth  In  the  ELF  band  (3  cps-3  kos) .  The 
basis  for  the  method  devised  was  to  use  worldwide  thunderstorm  activity  maps 
from  the  Geophysics  Handbook  for  Air  Force  Designers,  the  diurnal  thunder 

4 

probability  distribution  which  peaks  about  1600  hours  local  time,  and  the 


2.  Galejs,  J.,R.V.,Row,  Propagation  of  ELF  Waves  Below  an  Inhomogeneous 
Anisotropic  Ionosphere,  Research  Report  334.  10  April  1963.  Applied  Re¬ 
search  Laboratory,  Syivanla  Electronic  Systems.  To  appear  In  IEEE  Trans¬ 
actions  Prof.  Tech.  Grp.  Antennas  and  Propagation,  Jan.  1964. 

3.  Galejs,  J.,  ELF  and  VLF  Waves  Below  an  Inhomogeneous  Anisotropic  lonos- 

fhere.  Research  Report  350,  9  Xugust  196^,  Applied  Research  Laboratory, 
ylvania  Electronic  Systems. 

4.  Williams,  J.C.,  "Thunderstorms  and  Very-Low  Frequency  Noise",  Ph.D. 

Thesis,  Division  of  Engineering  and  Applied  Physics,  Harvard  University, 
June  1959. 
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spectrum  of  the  average  return  lightning  stroke^,  together  with  Galejs®  model 
of  propagation  below  an  exponentially  tapered  conductivity  versus  height 
profile  of  the  ionosphere,  to  predict  the  RMS  noise  spectrum  of  the  vertical 
electric  field.  Because  there  was  Interest  in  frequencies  well  above  the 
resonance  region  (5-35  cps)  it  was  necessary  to  extend  the  individual  dipole 
moment  stroke  spectrum  out  to  at  least  1  kc.  It  was  discovered  in  the  course 
of  attempting  this  that  the  available  data  on  stroke  current  waveforms  and 
statistics  is  insufficient.  Consequently  the  stroke  spectrum  was  taken  arbi¬ 
trarily  to  be  flat  out  to  1  kc.  Data  from  Ushuala,  Argentina  and  the  Panama 
Canal  Zone  (averaged  over  several  days)  were  used  to  normalize  the  absolute 
level  of  the  noise  at  a  frequency  of  35  cps.  Figure  5  shows  a  plot  of  the 
predicted  spectrum  at  Ushuala  and  Figure  6  the  predicted  spectrum  at  Malta. 
Originally  it  was  intended  to  have  DECO  measure  noise  simultaneously  at 
three  widely  separated  sites  and  compare  the  predictions  from  one  to  the 
other.  Such  measurements  were  made  in  June  1963  at  Malta,  Boulder,  Colorado 
and  Fairbanks,  Alaska.  However,  the  data  has  not  been  processed  up  to  this 
time. 

It  was  recognized  that  sources  near  an  observer  or  his  antipode  are  in¬ 
correctly  weighted  by  the  usual  simple  asymptotic  expansion  of  the  function 
P^(-cos  0) ,  which  was  used  in  the  above  computations.  Work  was  initiated  to 
develop  a  more  satisfactory  approximation  in  these  regions.  Lack  of  funds 
have  not  permitted  the  reworking  of  the  numerical  computations  including  this 
correction.  A  rough  calculation  has  been  made  where  the  contribution  of  noise 
sources  (lightning  strokes)  within  50  degrees  of  the  observer  and  his  antipode 
have  been  neglected.  Figure  7  shows  a  plot  of  the  calculated  diurnal  varia¬ 
tion  of  the  noise  at  frequencies  of  8,  14  and  20  cps  for  an  observer  at 
Ipswich,  Massachusetts.  Figure  8  shows  the  corresponding  experimental  ob¬ 
servations  of  Balser  and  Wagner.  In  view  of  the  crudeness  of  the  theoretical 
model  the  comparison  is  encouraging. 

The  thunderstorm  activity  maps  used  in  this  task  are  compiled  from  a 
variety  of  sources  and  refer  only  to  thunder  heard.  It  is  not  clear  a  priori 


5.  Raemer,  H.R.,  "On  the  Spectrum  of  Terrestrial  Radio  Noise  at  Extremely- 
Low-Frequencies , "  Jour.  Res.  NBS  65D,  581,  Nov.  1961. 

6.  Galejs,  J.,  "Terrestrial  Extremely-Low-Frequency  Noise  Spectrxim  in  the 
Presence  of  Exponential  Ionospheric  Conductivity  Profiles,"  Jour.  Geophys. 
Res.  66,  2787,  Sept.  1961. 
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that  there  is  a  one-to-one  correspondence  between  the  surface  density  of  the 
average  vertical  dipole  moment  strength  of  the  initiating  lightning  strokes 
and  thunder  Itself.  If  indeed  lightning  is  the  dominant  source  of  ELF  noise 
in  the  cavity  resonance  frequency  region  then  it  should  be  possible  to  work 
the  prediction  technique  in  reverse,  so  to  speak,  and  to  construct  maps  of 
lightning  dipole  density.  Work  was  started  on  determining  a  suitable  analyt¬ 
ical  -  computational  technique  for  processing  measurements  from  a  grid  of 
observing  stations  which  measure  the  RMS  or  average  noise  at  the  same  fre¬ 
quencies  simultaneously.  This  work  is  necessary  to  determine  how  many  sta¬ 
tions  and  frequencies  would  be  necessary  to  achieve  a  given  mapping  coverage 
of  the  world.  If  such  a  scheme  were  feasible  it  could  result  In  a  much 
better  picture  of  worldwide  thunderstorm  activity  and  serve  as  the  basis  of  a 
better  noise  prediction  service  at  VLF  and  possible  HP.  Unfortunately  no 
conclusive  results  were  available  at  the  time  of  writing. 

GUIDED  WAVES  NEAR  THE  AIR-IONOSPHERE  BOUNDARY 

A  Study  was  undertaken,  by  a  consultant,  into  the  possibility  of 
launching  guided  ELF  or  VLF  waves  from  a  source  near  the  air-ionosphere 
boundary.  This  would  pertain  to  a  transmitter  in  a  low  orbit  satellite,  to 

upward  lightning  strokes  and  possible  to  the  propagation  of  micropulsations. 

7  8 

The  work  is  reported  in  two  technical  reports.  ’ 


7.  Seshadrl ,  S.R.,  A.  Hessel,  Radiation  from  a  Source  Near  a  Plane  Interface 
Between  an  Isotropic  and  a  Syrotroplc  faielectrlc.  Research  Report  339. 
Applied  Research  Laboratory,  Sylvanl a  Electronic  Systems,  27  May  1963. 

8.  Seshadrl,  S.R.,  Radiation  from  Electromagnetic  Sources  in  a  Plasma. 
Research  Report  344 ,  Applied  Research  taooratory ,  Sylvanla  Electronic 
Systems,  5  June  1963. 


5 


F-0154-2 


REPORTS  AND  PUBLISHED  PAPERS 
ENGINEERING  NOTE  163,  J.  Galejs 

"Resonator  Type  Oscillations  Between  the  Earth  and  Ionosphere  and  Effects 
of  Ionosphere  Height  and  Conductivity  Changes  Between  Day  and  Night 
Hemispheres,"  5  July  1960. 

APPLIED  RESEARCH  MEMORANDUM  219,  J.  Galejs 

"Propagation  of  VLF  and  ELF  Radio  Waves  Excited  by  a  Horizontal  Magnetic 
Dipole,"  16  August  1960. 

APPLIED  RESEARCH  MEMORANDUM  223  Part  I  and  II,  J.  Galejs 
"Low  Frequency  Slot  Antennas,"  9  September  1960. 

ENGINEERING  NOTE  181,  H.  Raemer 

"The  Effect  of  Induced  Polarization  on  Propagation,"  28  September  1960. 

APPLIED  RESEARCH  MEMORANDUM  234  Part  I,  H.  Raemer 

"On  the  Spectrum  of  Terrestrial  Radio  Noise  at  Extremely  Low  Frequency," 

7  February  1961. 

APPLIED  RESEARCH  MEMORANDUM  240,  j.  Galejs 

"ELF  Waves  In  the  Presence  of  Exponential  Ionospheric  Conductivity 
Profiles,"  11  April  1961. 

ENGINEERING  NOTE  221,  J.  Galejs 

"Terrestrial  Extremely-Low-Frequency  Noise  Spectrum  In  the  Presence  of 
Exponential  Ionospheric  Conductivity  Profiles,"  25  April  1961. 

ENGINEERING  NOTE  234,  R.  Row 

"On  the  Resonant  Frequencies  of  the  Schmann  Modes , "  30  June  1961 . 
ENGINEERING  NOTE  243,  R.  Row 

"Tensor  Conductivity  and  ELF/VLF  Indices  of  Refraction  for  the  D-Region," 
31  July  1961. 

ENGINEERING  NOTE  280,  R.  Row 

"Vector  Wave  Equations  for  E  and  H  In  Inhomogeneous  Anisotropic  Media," 

6  December  1961 . 

ENGINEERING  NOTE  281,  J.  Galejs 

"Further  Note  of  Terrestrial  ELF  Propagation  In  the  Presence  of  an 
Isotropic  Ionosphere  with  an  Exponential  Conductivity  Height  Profile," 

7  December  1961. 

ENGINEERING  NOTE  294,  D.  Winter 

"The  Solution  of  a  Two-Dimensional  Analogue  of  the  Sommerfeld  Dipole 
Problem,"  9  January  1962. 


6 


F-0154-2 


ENGINEERING  NOTE  298,  D.  Winter 

"The  Magnetic  Field  In  the  Sea  Due  to  a  Two-Dimensional  Dipole  Array  at 
the  Surface,"  22  January  1962. 

ENGINEERING  NOTE  335,  J.  Galejs 

"Transverse  Propagation  In  an  Exponentially  Varying  Anisotropic  Medium," 

1  June  1962. 

RESEiWCH  REPORT  329,  D.  Winter 

"Low  Frequency  Radio  Propagation  Into  a  Moderately  Rough  Sea," 

5  March  1963. 

RESEARCH  REPORT  334,  J.  Galejs  and  R.  Row 

"Propagation  of  ELF  Waves  Below  an  Inhomogeneous  Anisotropic  Ionosphere," 
10  April  1963. 

RESEARCH  REPORT  350,  J.  Galejs 

"ELF  and  VLF  Waves  Below  an  Inhomogeneous  Anisotropic  Ionosphere," 

9  August  1963. 


Work  partially  supported  under  Nonr  3185(00) 

Research  Report  339,  S.  R.  Seshadrl,  A.  Hessel 

"Radiation  from  a  Source  Near  a  Plane  Interfact  Between  an  Isotropic 
and  a  Gyrotroplc  Dielectric,"  27  May  1963. 

RESEARCH  REPORT  344,  S.  R.  Seshadrl 

"Radiation  from  Electromagnetic  Sources  In  a  Plasma,"  5  June  1963. 


PUBLISHED  PAPERS 

1.  "Excitation  of  VLF  and  ELF  Radio  Waves  by  a  Horizontal  Magnetic  Dipole," 

J.  Galejs,  Jour.  Res.  NBS  65D,  305,  May-June  1961. 

2.  "Effect  of  Underground  Induced  Polarization  of  ELF  Propagation," 

H.  Raemer,  Jour.  Geophys.  Res.  66.  1596,  May  1961. 

3.  "On  the  Extremely  Low  Frequency  Spectrum  of  Earth-Ionosphere  Cavity 
Response  to  Electrical  Storms,"  H.  R.  Raemer,  Jour.  Geophys.  Res.  66, 

1580,  May  1961. 

4.  "ELF  Waves  In  the  Presence  of  Exponential  Ionospheric  Conductivity  Pro¬ 
files,"  J.  Galejs,  IRE  Trans.  Ant.  Prop.  AP-9 .  554,  Noveitdser  1961. 

5.  "On  the  Spectrum  of  Terrestrial  Radio  Noise  at  Extremely  Low  Frequencies," 
H.  R.  Raemer,  Jour.  Res.  NBS.  65D.  581,  November -December  1961. 

6.  "Terrestrial  Extremely  Low  Frequency  Noise  Spectrum  In  the  Presence  of 
Exponential  Ionospheric  Conductivity  Profiles,"  J.  Galejs,  Jour.  Geophys. 
Res . .  66 .  2787,  September  1961. 


7 


F-0154-2 


7.  "A  Further  Note  on  Terrestrial  Extremely  Low  Frequency  Propagation  in  the 
Presence  of  an  Isotropic  Ionosphere  with  an  Exponential  Conductivity 
Height  Profile,"  J.  Galejs,  Jour.  Geophys.  Res.  67,  2715,  July  1962. 

8.  "Terrestrial  Extremely  Low  Frequency  Propagation  in  the  Presence  of  an 
Isotropic  Ionosphere  with  an  Exponential  Conductivity  Height  Profile," 

J.  Galejs,  Proceedings  of  the  International  Conference  on  the  Ionosphere, 
London,  July  1962. 

9.  "On  the  Electromagnetic  Resonant  Frequencies  of  the  Earth  Ionosphere 
Cavity,"  R.  Row,  IRE  Trans.  Ant.  Prop.  AP-IO,  766,  November  1962. 

10.  "Low  Frequency  Radio  Propagation  into  a  Moderately  Rough  Sea,"  D.  F. 
Winter,  Jour.  Res.  NBS.  67D,  551,  September -October  1963. 


8 


F-0154-2 


3-1-0372 


Figure  3.  Normalized  Propagation  Constants  -  Nighttime  -  With  Ions  Considered 
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Figure  5.  ELF  RMS  Noite  Field  -  Ushuoia,  Argentine  -  1500  GMT 
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Figure  6.  ELF  RMS  Noitc  Field  -  Itlartd  of  Malta  -  1500  GMT 
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Figure  7.  Cglculoted  Diurnal  Variation  of  Noiie  Pbwer  at  the  Peak  of  the  Resonance  Modes  - 
Wintertime  -  Boston,  Massachusetts 
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UNIVERSAL  TIME 

Figure  8.  MMtured  Diurnal  Variation  of  Noiit  Pbwar  at  tha  Poak  of  tha  Rotonanco 
Modot  -  Ipiwieh,  MosiochuMttt  -  Fabruary,  1961 
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